Quantum dynamics of two bosons in an anharmonic trap: Collective vs internal excitations 
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This work deals with the effects of an anharmonic trap on an interacting two-boson system in one dimension. 
Our primary focus is on the role of the induced couphng between the center of mass and the relative motion as 
both anharmonicity and the (repulsive) interaction strength are varied. The ground state reveals a strong local- 
ization in the relative coordinate, counteracting the tendency to fragment for stronger repulsion. To explore the 
quantum dynamics, we study the system's response upon (i) exciting the harmonic ground state by continuously 
switching on an additional anharmonicity, and (ii) displacing the center of mass, this way triggering collective 
oscillations. The interplay between collective and internal dynamics materializes in the collapse of oscillations, 
which are explained in terms of few-mode models. 
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I. INTRODUCTION 

Ultracold bosonic atoms have become a popular tool for 
simulating and understanding fundamental quantum phenom- 
ena (see Refs. 1 1, 2, 3, 4] for an overview). The key reason is 
that, drawing on their interaction with electromagnetic fields, 
cold atoms permit a tremendous degree of tunability of both 
external (trapping) forces and atomic interactions. Typically, 
however, studies of these kinds of systems rest on two essen- 
tial premises: harmonicity of the trap, and mean-field interac- 
tions. The first assumption is in fact very natural and easily 
fulfilled experimentally, at least near a well-defined minimum 
of the potential. In fact, the separability of the center-of-mass 
motion for identical particles facilitates many theoretical in- 
vestigations and has important ramifications for the dynam- 
ics, such as the existence of undamped dipole oscillation. The 
obvious question as to how additional anharmonic forces in- 
fluence the eigenstates as well as the dynamics has attracted 
little attention. In ref. [5], e.g., the effect of anharmonic terms 
on the spectrum of two unequal atoms is studied, whereas 
Ref. [6] addresses the collapse of collective oscillations, but 
from a mean-field perspective. 

Within that mean-field approximation, all atoms are as- 
sumed to coherently populate a single-particle orbital, which 
is governed by the Gross-Pitaevskii equation While this 
is justified in the limit of atom numbers 3> 1 and suf- 
ficiently weak interactions, it neglects two-body effects that 
become relevant for few-body systems, and in particular in 
the limit of strong transversal confinement, when the sys- 
tem becomes effectively one-dimensional (ID). In that case, 
it has been shown ^ that the effective coupling strength can 
be tuned at will, making it possible to change from a non- 
interacting to a strongly repulsive system. In that highly cor- 
related borderline case, the bosons are said to fermionize [Sj] 
as they acquire partly counterintuitive properties similar to 



an ideal gas of fermions. This, among other rationales, has 
sparked many efforts toward understanding few-boson sys- 
tems Hill [Bill]. 

The subject of this paper is to carry out a systematic ab ini- 
tio analysis of anharmonic effects on the ground state as well 
as on the quantum dynamics of interacting bosons. We focus 
on the fundamental two-atom case, which already reveals a 
rich variety of features, thus paving the way to systems with 
more atoms. It is solved utilizing an exact and flexible multi- 
configurational approach, which has proven very fruitful in 
previous studies on the fermionization of ensembles with up 
to SIX atoms Here the ground state is the nat- 

ural starting point for an investigation into the nature of the 
coupling between the center of mass and the relative coordi- 
nate. Equipped with this insight, we probe the impact of that 
coupling on the dynamics via two general schemes: First, we 
switch on an additional anharmonic trapping force, this way 
exciting the (uncoupled) harmonic ground state. Second, we 
excite the collective motion by elongating the center of mass 
and follow the induced dynamics in the relative motion. 

This paper is organized as follows: After a brief introduc- 
tion of the theoretical model in Section|II] we present our com- 
putational method in Sec.|III] The subsequent sectionHVldeals 
with the ground state of the two-boson system in an anhar- 
monic trap. In Sec. |V] we study the dynamics after excitation. 



II. MODEL 

In this work we investigate a system of two interacting 
bosons under external confinement. These particles, repre- 
senting atoms, are taken to be one-dimensional (ID). More 
precisely, after integrating out the transverse degrees of free- 
dom and rescaling to dimensionless units, we aiTive at the 
model Hamiltonian (see 1,14,1 for details) 



H = hi + h2 + V{xi 



(1) 
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where hi — ^pf + U{xi) is the one-particle Hamiltonian with 
a trapping potential U, while V is the effective two-particle 
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interaction V{x) = g5cr{x). Here we concentrate on re- 
pulsive foces, g > 0- The well-known numerical difficulties 
due to the spurious short-range behavior of the standard delta- 
function potential S{x) are alleviated by mollifying it with the 
normahzed Gaussian 
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which tends to 6{x) as ct ^ in the distribution sense. We 
choose a fixed value ct = 0.1 as a trade-off between smooth- 
ness and a range much smaller than the average particle dis- 
tance. 

We assume a harmonic trap superimposed by an anhar- 
monic potential with linear and quartic terms: U{xi) = 
+ KXi + \x^. The linear part of the potential causes a 
displacement of the trap center (controlled by a displacement 
factor k), while the quartic terms squeeze the trap (controlled 
by an anharmonicity factor A). By time-dependent variation 
of both parameters later on in this work, the trap can be dis- 
torted so as to cause excitations in the two-boson system. 

For studying the effects of the coupling between the cen- 
ter of mass and the internal relative dynamics, we perform a 
coordinate transformation to the center-of-mass frame of ref- 
erence: R — {xi + X2) and X — xi — X2- Finally, we arrive 
at the following Hamiltonian 
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After the transformation to the center-of-mass frame, the in- 
teraction depends only on the relative coordinate x, whereas 
the shift of the trap only affects the center of mass coordinate 
R. Moreover, for a quartic anharmonicity (A > 0) an Rx- 
coupling term appears in the Hamiltonian. 



III. COMPUTATIONAL METHOD 

Our goal is to investigate the ground state and the dynamics 
of the system introduced in Sec. |II] for all relevant interac- 
tion strengths in a numerically exact, i.e., controllable fash- 
ion. Our approach relies on the Multi-Configuration Time- 
Dependent Hartree (MCTDH) method [17], primarily a wave- 
packet dynamics tool known for its outstanding efficiency in 
high-dimensional applications. To be self-contained, we will 
provide a concise introduction to this method and how it can 
be adapted to our purposes. 

The underlying idea of MCTDH is to solve the time- 
dependent Schrodinger equation 



(4) 



as an initial-value problem by expansion in terms of direct (or 



Hartree) products $,/ = (ysj^-* 
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*(.,i)=^Aj(i)<i>j(.,i), 



(5) 



using a convenient multi-index notation for the configurations, 
J — {ji, . . . ,jf), where / denotes the number of degrees of 

freedom. The (unknown) single-particle functions ^p^^^ (Jk = 
1, . . . , 71k) are in turn represented in a fixed primitive basis 
implemented on a grid. In our work, we consider systems with 
two degrees of freedom and use rifi = Ux = orbitals. The 
grid spacing, in our case, should of course be small enough 
for the relative coordinate to sample the interaction potential, 
whereas the center-of-mass coordinate is not that sensitive. In 
this light, we employ a discrete variable representation with 
95 and 165 grid points for the center-of-mass coordinate and 
the relative coordinate, respectively. 

Note that in the above expansion, not only the coefficients 
Aj are time dependent, but so are the Hartree products 
Using the Dirac-Frenkel variational principle, one can derive 
equations of motion for both Aj,(pj |fl7j]. Integrating this 
differential-equation system allows one to obtain the time evo- 
lution of the system via (|5]l. Let us emphasize that the concep- 
tual complication above offers an enormous advantage: the 
basis {$j(-, f)} is variationally optimal at each time t. Thus 
it can be kept fairly small, rendering the procedure very effi- 
cient. 

The MCTDH approach ifTsll . which we use, incorporates a 
significant extension to the basic concept outlined so far The 
so-called relaxation method [19] provides a way to not only 
propagate a wave packet, but also to obtain the lowest eigen- 
states of the system, ^fQ. The key idea is to propagate some 
wave function ^I*" by the non-unitary e~^^ (propagation in 
imaginary time.) As t — > 00, this exponentially damps out 
any contribution but that stemming from the true ground state 
like e^^"^. In practice, one relies on a more sophisticated 
scheme termed improved relaxation IqBl I21I1 . which is much 
more viable especially for excitations. Here ('I'|i/|'I') is mini- 
mized with respect to both the coefficients Aj and the orbitals 
(fij . This leads to (i) a self-consistent eigenvalue problem for 
{{^j\H\^k)), which yields Aj as 'eigenvectors' , and (ii) 
equations of motion for the orbitals (pj based on certain mean- 
field Hamiltonians. These are solved iteratively by first diago- 
nalizing for Aj withfixed orbitals and then 'optimizing' ipj by 
propagating them in imaginary time over a short period. That 
cycle will then be repeated until convergence is achieved. 



IV. GROUND STATE OE TWO-BOSON SYSTEMS IN 
ANHARMONIC TRAPS 

Before studying the quantum dynamics of the interacting 
two-boson system in an anharmonic trap, we examine its 
ground state in dependence of the interaction strength g > 
and the anharmonicity A. By increasing A, and thus squeez- 
ing the trap, we can force the particles to reduce their distance, 
this way counteracting the tendency to fragment as the interac- 
tion strength g ^ 00. Figure [T] depicts one-particle densities 
pi{x) — J dR I (i?, x) p of the relative coordinate for differ- 
ent interaction strengths g and anharmonicities A in an undis- 
placed (k = 0) anharmonic trap. Obvioulsy [cf. Fig. [n^)]. 
for non-interacting particles, the probability density is peaked 
at a; = 0, i.e., both atoms remain in the trap center regardless 
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of whether there is another one already. For increasing repul- 
sion [Figs. [Ilb-d)], the mean distance between the particles 
grows, and more and more a pronounced fragmentation in the 
density occurs, indicating the inhibition of both particles be- 
ing located at the same point in space. At g — 20.0 (cf. fig. 
[T](d)), the fragmentation is almost complete. By contrast, the 
squeezing of the trap caused by the anharmonic term A > 
leads to a reduction of the fragmentation in the density. Both 
effects make for a stronger localization of the two particles. 



(a) 



(b) 
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Furthermore, to illuminate the role of the coupling between 
the center of mass (i?) and the relative coordinate {x), we con- 
sider the same situation as before, but we artificially switch off 
the coupling term 3XR^x^ in the Hamiltonian in Eq. (|3]l. For 
a special trap with anharmonicity A = 6.0 the corresponding 
two-particle densities p2{R,x) = a;)p with and with- 

out i?x-coupling are shown in Fig. |2] for different interaction 
strengths (note that only for two degrees of freedom, as in our 
case, the two-particle density is equal to the absolute square of 
the wave function). Comparing the true results [Figs. |2|a-c)] 
with those where the i?a;-coupling term has been artificially 
excluded [Figs. |2|d-f)], it becomes clear just how that cou- 
pling suppresses configurations in which both the center of 
mass and the relative coordinate are large. In other words, 
it is even less likely for the atoms be to be far from the trap 
center and simultaneously far from each other. This is vis- 
ible in the two-particle densities as an unphysical flattening 
in the affected regions in the absence of the Rx term [Figs. 
|2d-f)], in contrast to the full result. The effect is more pro- 
nounced for larger anharmonicities A, when the coupling term 
becomes more relevant (although it is still dominated by the 
quartic terms like R'^ especially for large values of R and x). 

To measure the mean relative distance between the two par- 
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Figure 1: Ground state of the two-boson system in an anharmonic 
trap with different anharmonicity factors A. Shown are the one- 
particle densities for the relative coordinate pi{x) for interaction 
strengths of (a) g = 0.0, (b) g = 1.0, (c) g = 6.0, (d) g = 20.0 



Figure 2: (a) - (c): Two-particle densities p2{R,x) of the ground 
state of the two-boson system in an anharmonic trap with anhar- 
monicity factor A = 6.0 and (a) g = 0.0, (b) g = 1.0, (c) g = 20.0. 
— (d) - (f) same, but without 7?a:-couphng term in the Hamiltonian. 
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Figure 3: (a) Mean relative distance of the two particles and (b) en- 
ergy of the groundstate of the two-boson-system both as a function 
of the anharmonicity factor A and the interaction strength g. 

tides, we look at the quantity (x^), which is equal to the 
square of the standard deviation Ax, as (x) — due to permu- 
tation symmetry. It turns out [Fig. Oa)] that with increasing 
interaction strength (for a constant anharmonicity of the trap 
potential) the mean relative distance first rises quickly and 
then becomes almost constant because of the limited range 
of the interaction. By contrast, a larger anharmonicity in the 
trap potential and thus a stronger confinement of the particles 
leads to a decrease in the (mean-square) relative distance, as 
expected. Figure [3lb) shows the energy of the ground state 
as a function of the interaction strength g, plotted for differ- 
ent anharmonicity factors A. For a constant anharmonicity the 
energy first increases with growing interaction strength, then 
more and more saturates because of the finite range of the in- 
teraction. For a larger anharmonicity, the energy of the ground 
state is higher because of the increase in potential energy. Ow- 
ing to the squeezing, it is harder for the atoms to move apart 
and isolate each other, this way shifting the fermionization 
limit to larger g. 



ing function Q{t) which was chosen as 

e(i) = i(l±tanh[e(t-r)]). (6) 

For all intents and purposes, this function is time independent 
except for a short while controlled by the parameters ^ (the in- 
verse duration of the switching) and r (responsible for the mo- 
ment in time of switching). The positive (negative) sign cor- 
responds to switching on (off) the potential terms. This way, 
the Hamiltonian is essentially time-dependent only during the 
switching process. Before that, the system is in its ground 
state; after the excitation (for t — r > 3 >• 1/0, it is propa- 
gated by the time-independent Hamiltonian with k = 0. 
If the duration of the switching is sufficiently long compared 
to the time scales of the system, the excitation process oc- 
curs adiabatically and the system goes over into to its mo- 
mentary ground state (for a system with parameters g = 20.0 
and A = 0.1 this is the case for ^ < 0.5). In the following we 
use the values ^ = 3.0 and t = 2.0. 

B. Excitation by switching on an anharmonic potential 

For the first excitation scheme, we prepare the two-boson 
system in the ground state of a purely harmonic trap (A ~ 0). 
Then, over a short while characterized by the parameters ^ 
and T (cf. eq. (|6]l), an additional quartic potential is switched 
on (with an asymptotic value A > 0). After that, the system — 
now again described by the time-independent Hamiltonian (O 
— is propagated for some time t, and its quantum dynamics is 
studied. 



V. QUANTUM DYNAMICS 

After having explored the ground state, we now want to 
consider excitations of the interacting two-boson system. Al- 
together, we apply two methods of excitation, both realized 
by a distortion of the external trap. In the first scheme, de- 
scribed in Sec. IV B| we switch on a quartic anharmonic po- 
tential in addition to a purely harmonic trap; in the other 
case (Sec. IV Cb . we displace the trap center by adding a lin- 
ear potential to a quartic anharmonic trap. In both cases we 
demonstrate that collapse and revival of oscillations in the 
relative coordinate occur if the excitation is performed non- 
adiabatically, that is, if the duration of the switching process 
is short compared to the relevant time scales of the system. 
Moreover, we analyze the excitation spectra in order to ex- 
plain the dynamics within a few-mode model. 



A. Switching process 

Let us first briefly discuss the switching procedure itself. 
Technically the switching is implemented by multiplying the 
potential terms to be switched on or off with a specific switch- 



1. Collapse and revival of oscillations 

The excitation process outlined above alters the shape of 
the trap, but without affecting the parity symmetries for both 
R and x; hence, (R) — for all times (recall that trivially 
(x) = by permutation symmetry alone). However, the 
squeezing excites breathing oscillations both in the center-of- 
mass width AR (omitted here) as well as in the internal mo- 
tion, (Ax)^ — (a;^), which experiences collapses and revivals 
[cf. Figs. Eta-d)]. There are time spans in which the inter- 
atomic distance is rapidly changing, while at other times the 
system is almost at rest or only oscillating very slowly. The 
time between two collapses of the dynamics depends on the 
excitation process, i.e., on the amount of energy added by the 
anharmonic potential as well as on the interaction strength be- 
tween the two particles. It turns out that the stronger the added 
anharmonic potential (cf. figs. |4](a/b),(c/d)) and the stronger 
the interaction (cf. figs. |4](a/c),(b/d)) are, the shorter is the 
time between two collapses. 

2. Excitation spectra and symmetry analysis of excited states 

For a better understanding of the effects of the excita- 
tion, it is natural to analyze the contributing excited states of 
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Figure 4: Dynamics of the two-boson system after switching on a 
quartic anharmonicity. Square of the mean relative distance of the 
two particles {x^) for (a) g = 1.0, A = 0.1, (b) g = 1.0, A = 0.5, 
(c) g = 20.0, A = 0.1, (d) g = 20.0, A = 0.5 and — (e) g = 20.0, 
A = 0.1 (without _Ra;-coupling term) 



the stationary Hamiltonian H{t ^ 1/0' i-^ ' including the 
full anharmonic potential and without the displacement term 
(A > 0, K — 0). We describe the system after excitation by 
expanding its wavefunction in terms of the eigenstates 



^I^, of that Hamiltonian, 



(7) 



where H'i>j{R, x) — ujj'^ j{R, x) and ujq denotes the ground 
state energy. The analysis of the excitation spectrum — ob- 
tained via Fourier transformation of the autocorrelation func- 
tion a{t) = (vl'(O)l^(i)) — provides us with the values for 
LOj . On the other hand, the exact eigenvalues can also be ob- 
tained by diagonalizing the Hamiltonian and then be related 
to the peaks in the spectrum a(w). 

This procedure yields that the ground state '^q gives the 
main contribution (cf. fig. |5]l . Furthermore, not all but only 
certain higher states are populated in For very large inter- 
action strengths {g > 10), apart from the ground state in the 
anharmonic trap, only the eigenstates j — 4, 6, 12, 14, 16 are 
excited. For smaller interaction strengths (g < 10), in turn, 
the lowest eigenstates j = 4, 5, 12, 13, 16 are excited. Gen- 
erally, higher states are less and less excited while the ground 
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Figure 5: Excitation spectra for (a) g — 20.0, A = 0.1, - (b) g = 
20.0, A = 1.0. The numbers indicate the excited states, denotes 
the ground state. 



state contains the main contribution. The more energy is put 
into the system by the distortion of the harmonic trap (cf. fig. 
|5](a) and (b) for g = 20.0 and A = 0.1 and 1.0, respectively) 
the more strongly energetically higher states are excited. The 
same is true for fast excitation processes, which likewise de- 
posit more energy into the system. For small distortions of the 
harmonic trap (A < 0.1) and an adequate duration of switch- 
ing {£, ^ 3.0), only two higher states (the 4*'' and 6*'' for 
large, the 4*'* and for small interaction strengths) are con- 
siderably excited while all other states have only a negligible 
amplitude. 

The fact that not all of the lowest eigenstates are excited 
hints at selection rules due symmetries of the system. Indeed, 
as is obvious from Eq. [3] (for k = 0), the Hamiltonian con- 
serves both the center-of-mass parity and the relative parity, 
so only states having the same R- and a; -parities as the initial 
state can be excited. A closer analysis confirms that only the 
states mentioned above have even i?-parity as well as even x- 
parity and thus obey the same symmetries as the initial state, 
i.e., the harmonic ground state. It should be mentioned that in 
our way of counting the excited states, we include all states 
having odd a; -parity (e. g. the P* or the 3'"'' excited state), 
which would represent a non-bosonic system. 



3. Few-mode model 

The dynamics of the system can be fitted very well to a 
simple few-mode ansatz. For this, we revert to the expansion 
(Q of the wavefunction, restricted to to + 1 eigenstates of the 
Hamiltonian. Solving the Schrodinger equation fixes the time- 
dependent coefficients in the most general form as Cj{t) = 
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Cj e~*"j *^"^j , where the arbitrary phases dj will be used in 
the fits to account for the duration of the switching process. 
The values for Cj G R can be obtained from the height of the 
corresponding peak in the spectrum. 

With this ansatz, expectation values like ^x^) can be easily 
calculated, yielding a constant, time-independent part and a 
dynamic part made up of a sum over harmonic oscillations: 

^ rn m 
j=0 3^k=0 

where the frequencies are given by ujj^k = i^j — ^k, the phases 
by dj^k = dj — dk and the amplitudes by 

aj,fc = 2cfeCj(«'fc|a;2|^'j). (9) 

In an analogous way, also expectation values describing the 
center-of-mass dynamics like (i?) or {E?) can be calculated. 
They share the frequencies of the modes uoj,k with the dynam- 
ics of the relative motion but the amplitudes aj,k are different. 

Within this few-mode model, we can now explain the dy- 
namics of the two-boson system after excitation via distor- 
tion of the external trap. The frequencies Wj fc of the modes 
contributing to the dynamics are determined by the frequen- 
cies of the contributing excited states. For weak excitations 
as mentioned above, only two excited states are consider- 
ably populated, and in this case it is expected that only three 
modes determine the dynamics of the system (for large inter- 
action strength these are the ones with frequencies wo,4, wo,6 
and ^4.6)- Even more so, since the amplitude of the third 
mode (with frequency 1JJ4.6) is much smaller than the ampli- 
tudes of the two other modes and has a much larger timescale 
(w4 6 ^ ijJo.4, ^vi.q), practically only two modes are expected 
to determine the dynamics of the two-particle system for weak 
excitations: oj^^^ = a;4 — cjq and cj'^^) = ujq — ujq. In this case 
we expect a simple beat behavior in the dynamics, where the 
beat period determining the time between two collapses is in- 
versely proportional to w'^' — w'^' = tjg ~ ^4 for large inter- 
action strengths (cjs — oj/^ for small ones). The condition for a 
complete collapse is the equality of the amplitudes of the two 
contributing modes (ao,4 ~ ao.e)^ which is almost fulfilled 
in the relative dynamics but not valid for the center-of-mass 
breathing oscillations. So, by the way of excitation, e.g. the 
extent of anharmonicity added to the potential, it is possible 
to control the internal dynamics of the system, e.g., the time 
between two collapses of the oscillations. 

Figure |6] (a) confirms that the presented few-mode model 
(including 2 modes) fits the simulated data very well. Apart 
from that, the optimal fitted values for the frequencies of the 
modes coincide perfectly with the theoretically predicted val- 
ues for Wo, 4 and ujo^. However, a closer look reveals that for 
regions in which the relative dynamics has a small amplitude 
[Fig. |6lb)], minor deviations emerge which can be fixed by 
including more modes in the theoretical model. Within an 
adequate 4-mode model (by adding the next strongest modes 
with frequencies ti>4 g and cjg.ie), satisfactory agreement with 
the simulated data can be established. For stronger excita- 
tions, even more and higher states are considerably excited. 
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Figure 6: Fits of {x^) for a system with parameters g = 20.0 and 
A = 0.1. The open pentagons represent the exact simulation, (a) 
2-mode fit for time interval t G [10, 40] — (b) and (c): more detailed 
2-mode- and 4-mode fits, (b) for time interval t G [20, 26], (c) for 
time interval t £ [33, 39]. 

and thus even more modes contribute significantly to the dy- 
namics of the system. Then, of course, a pure 2-mode model 
no longer provides a good description for the dynamics of the 
system and for a satisfactory agreement between theory and 
simulation more modes must be included. 



4. Excitation without Rx-coupling term 

If the system is excited with an anharmonic potential in 
which artificially the i?a;-coupling term is eliminated, it turns 
out that the same states are excited (i.e. the new ground state, 
the 4*'\ 6*'* state and so on for larger interaction strengths) but 
some of the matrix elements emerging in equation (|9]) vanish 
(e.g. ao,4 = 0) so that in the end fewer modes contribute to 
the dynamics of the system. 

For weak excitations, the dynamics is dominated by a sin- 
gle mode with frequency wo,6 [cf. Fig. Hfe)]; thus there is 
no beat behavior anymore. Hence, in contrast to the full ex- 
citation (including the i?x-coupling term in the Hamiltonian) 
there is no collapse in the relative dynamics. This gives a strik- 
ing illustration of the role played by the coupling between the 
center-of-mass and the relative motion. 



C. Excitation by elongating tlie center of mass 

The second method of excitation examined in this paper 
uses the displacement of the trap center by a linear poten- 



7 





15 



20 



25 
t 



30 



35 



Figure 7: Dynamics in thie anhiarmonic trap after elongating thie cen- 
ter of mass. The open pentagons represent the exact simulation. 

(a) 3-mode fit of the center-of-mass position (R) for g = 20.0, 
A = 0.1, K = 1.0 

(b) and (c): Fits of the mean-square distance of the two particles (a::^) 
for g = 20.0, A = 0.1, k = 1.0 in the time interval t G [15, 35]: (b) 
3-mode fit, (c) 7-mode fit. 



tial. We start with the ground state of the two-boson system 
in an anharmonic trap whose center is displaced (Eq. ([3]) with 
A > 0, K > 0) and then continuously turn off the linear poten- 
tial (k 0) to revoke the displacement. With this procedure 
we excite collective oscillations: In fact. Fig. |3a) reveals how 
the center of mass oscillates about the trap center Moreover, 
since the anharmonic force couples the center of mass to the 
relative coordinate, also breathing oscillations in the relative 
motion are induced. These are displayed in the time evolu- 
tion of the widths Ax^ — (x^) in Fig. |7Jb). Again, we ob- 
serve collapses and revivals of the oscillations in the relative 
dynamics, whereas in the center-of-mass motion the ampli- 
tude of the oscillations periodically changes with time but the 
oscillations never die out completely. As both excitation tech- 
niques are similar, we can adopt much of the argumentation 
from the last subsection and focus mainly on the differences. 
Altogether, we find that the dynamics is more complicated 
than in the previous scheme. The reason for that will be laid 
out in the following two paragraphs. 



1. Excitation spectra and symmetry analysis 

An analysis of the spectra shows that, on top of the states al- 
ready excited with the previous excitation scheme, now more 
states contribute, the lowest ones being the 2"'', 8*'' and 10*'' 
state. However, along the same lines as before, this is readily 
explained in terms of the displacement term kR in the Hamil- 
tonian ([3j, which breaks the i?-parity (so that only the trivial 
x-parity is conserved during the excitation process). Hence, 



the 2 



nd oth 



and 10 state — all having odd i?-parity — are 



now also available. 



2. Few-mode model 

Equipped with this knowledge, the simulation can again be 
explained within the simple few-mode model introduced in 
Sec. IV B 3] The only thing is that, as now more states are ex- 
cited, also more modes significantly contribute to the dynam- 
ics, which is therefore more complicated as compared to the 
other excitation method. Especially the 2'"^ state is, next to 
the ground state, the most populated state after excitation, but 
also the 8"* and lO*'' state are more populated than e.g. the 
16*'*. As an example, for an interaction strength of g = 20.0, 
an anharmonicity of A = 0.1 and a displacement of k = 1.0, 
Figs. |7](b) and (c) illustrate that for the relative dynamics at 
least seven modes are necessary to have a satisfactory agree- 
ment of the simulation with the theoretical model. On the 
other hand, the dynamics of the center of mass can already 
be very well described in terms of only three modes having 
frequencies 0^0,2, 1^2,4 and see Fig. |2ta). 



VI. CONCLUSION AND OUTLOOK 

In this work we have investigated a repulsively interact- 
ing two-boson system in a one-dimensional anharmonic trap. 
We have investigated both the ground state of the system as 
well as its quantum dynamics upon excitation via distortion 
of the trap, with an eye toward the impact of the coupling 
between center of mass and relative motion. Our calculations 
are based on the numerically exact Multi-Configuration Time- 
Dependent Hartree method. 

Two different methods of excitation have been applied. In 
the first one we continuously switch on an anharmonic poten- 
tial in addition to an otherwise harmonic trap. This process 
induces breathing oscillations in the relative motion, which 
experience collapses and revivals. The dynamics can be ex- 
plained within a simple few-mode model based on the con- 
tributing excited states. For weak excitations, only two modes 
determine the dynamics of the system, while for stronger ex- 
citations more modes must be included. Just which states are 
excited depends on the symmetries of the initial state and thus 
on the method of excitation. With this knowledge it is possible 
to control the internal motion of the system — viz., the ampli- 
tude of the oscillations or the time between two collapses — by 
the external excitation via an adequate distortion of the trap- 
ping potential. We have also illuminated the role played by 
the coupling between the center of mass and the relative mo- 
tion by artificially excluding the corresponding terms in the 
Hamiltonian. Then, at least in the limit of weak excitations, 
one single mode dominates the oscillations, which are thus 
etirely undamped. 

In the second excitation scheme, we trigger collective os- 
cillations by displacing the center of mass, whose coupling 
to the relative coordinate leads to similar internal excitations 
as before. However, since the initial displacement breaks the 
center-of-mass parity, the dynamics is now more complex. 
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An obvious extension of our work would be the study of 
the dynamics of systems with more than two bosons, this way 
gaining insight into anharmonicity effects on a many-body 
level. One could also conceive more sophisticated excitation 
schemes so as to get detailed probes of the internal excitation. 
This may help not only understand experimental effects be- 
yond harmonic confinement, but also to actively control the 
interatomic dynamics. 
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